Anisotropic in-plane thermal conductivity of single-crystal YE^C^Os 
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We report measurements of the in-plane thermal conductivity (k) of YBa2Cu4C>8 (Y-124) single 
crystals in the temperature range 4K < T < 300A" and compare them with previous results on 
YBa2Cu306+a; (Y-123). For transport perpendicular to the CuO chains, K a (300K) ~ lOW/mK, 
and along the chains, Kt(300K) ~ 40W/mK. The temperature dependence of k for both transport 
directions is much stronger than in YBa2Cu306+a; (Y-123), indicative of substantially superior 
lattice heat conduction in Y-124, and resulting in maximum values for k at T ~ 20K exceeding 
200W/mK. The data imply a surprisingly large anisotropy in the lattice conduction. K a and Kb are 
enhanced in the superconducting state as in other cuprates. The magnitude and anisotropy of the 
enhancement are discussed and compared to those of Y-123. 

PACS numbers: 74.72.-h, 74.25.Fy, 74.25.Kc, 74.25.Dw 



The in-plane thermal conductivity (K a b) of cuprate 
superconductors has been studied extensively over the 
past several years. EJ Unresolved issues remain concern- 
ing the normal-state temperature dependence and the 
superconducting-state enhancement. Because of the low 
carrier density, heat conduction by the lattice accounts 
for more than half of the measured normal-state K a b in 
these materials. One would expect K a b(T) in high-quality 
crystals to reflect a dominant lattice contribution simi- 
lar to that of crystalline insulators, i.e. K a i, ~ 1/T at 
high-T, rising sharply to a maximum at low T.|-jThis be- 
havior is observed in Nd2- a; Ce a: Cu04(Nd-214)p but for 
all other cuprates investigated K a b is very weakly T de- 
pendent for T > T c . Whereas for Bi2Sr2CaCu0.8 struc- 
tural disorder (e.g., the Bi-0 layer modulation!]) imight 
explain its nearly glass-like lattice contribution,^ this 
would not seem applicable to YBa2Cu30g+a: (Y-123) and 
La 2 _a;Sr a ;Cu04 (La-214). Furthermore, K a b(T) for un= 
doped, insulating Y-123 and La-214 is unconventional 
suggesting that the weak T-dependence of K a b is generic 
to the latter materials and not directly related to 
the presence of free charge. An unidentified phonon 
damping mechanism appears to underlie_this behavior, 
possibly due to tocal lattice distortionsQ and/ice, mag- 
netic excitations B Regarding the enhancement pa ther- 
mal Hall conductivity experiments!!! on optimally-doped 
Y-123 imply that the phenomenon is largely electronic in 
origin, but less is known about underdoped compounds. 
The magnitude of the enhancement correlates with the 
superconducting pair densityBEj in Y-123 throughout 
the underdoped regime, and it is of interest to explore 
the generality of this behavior through studies of other 
cuprates. 

Here we report measurements of the in-plane thermal 
conductivity of YBa2Cu40s (Y-124). This material is of 
particular interest with regard to the issues mentioned 
above because of its structural similarity to Y-123, ab- 
sence of oxygen vacancies on the CuO chains, and the 



naturally underdoped state of the Cu02 planes. We find 
that the magnitude of n a (transverse to the chains) at 
T=300K is comparable to that of Y-123, but Kb is 3- 
4 times larger. A very large in-plane anisotropy in the 
lattice conduction is implied. Both K a (T) and Kb(T) be- 
have as in conventional crystalline insulators (like Nd- 
214); their maximum values (at T ~ 20K) are the highest 
reported for any cuprate and exceed those of Y-123 by an 
order of magnitude. These results imply a strong damp- 
ing of phonons by static or dynamic structural distortions 
as the source for the much weaker T dependence of K a b 
in Y-123. The superconducting-state enhancement for 
Y-124 is comparable in magnitude to that of T C =60K Y- 
123, consistent with the underdoped state of the planes in 
the double-chain compound. We discuss the anisotropy 
of the enhancement and its implications for the Lorenz 
number. 

Thermal conductivity measurements were performed 
on three single crystals grown Jay a high-pressure flux 
method as described previously.til Two of these (189 and 
259) were grown in AI2O3 crucibles yielding a slight Al 
contamination, and T c = 72K. Recent analyses indicateO 
that Al substitutes for 1-2% of the Cu(2) atoms in the 
Cu0 2 planes. As we discuss below, this light doping has a 
substantial effect on the thermal conductivity. The third 
crystal (315), grown in Y2O3 and without Al contami- 
nation, had T c — 80K. Typical crystal dimensions were 
0.8 x 0.5 x 0.05mm 3 , with the shortest dimension along 
the crystalline c axis. The a- and 6-axis electrical resis- 
tivities of similarly-prepared crystals have been discussed 
extensively elsewhere;^ typical values at T =300K are 
p a = 400/ificm and pi, — 130/iSlcm. The steady-state 
thermal conductivity measurements employed a fine- wire 
differential chromel-constantan thermocouple and minia- 
ture chip resistor as heater, both glued to the specimen 
with varnish or epoxy. The absolute accuracy of the mea- 
surements is ±20% due to uncertainty in the placement 
of the thermocouple junctions. No corrections for heat 
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losses (via radiation and conduction through the leads) 
have been applied; these are estimated to be ~ 10% 
near room temperature and < 2% at T < 120K. For 
specimens 259 and 315 gold contacts were vapor de- 
posited and their a-axis thermoelectric powers measured, 
yielding S(290K) = 36/iV/K and 16/xV/K, respectively. 
These values providp-cstimatcs of the hole concentration 
per planar Cu atom,Eilp ~ 0.09 and 0.11 for the Al-doped 
and undoped specimens, respectively. 

The results for n a and Kb are shown in Fig. [l]. 
The anisotropy in the normal state is substantial, with 
Kb/K a ~ 3—4 (Fig. ||). The normal-state T dependence in 
both crystallographic directions is close to 1/T for spec- 
imen 315 and is significantly weaker (~ 1/T 1 / 2 ) for the 
crystals containing Al. In all cases the T-dependence is 
substantially stronger than that of both T c = 92K and 
60K Y-123 (solid lines, Fig. |). At low T the character- 
istic dielectric maxima are observed near 20K. A change 
in slope is evident at T c [Fig. || (a)] for all specimens, 
though not appearing as prominent as in Y-123 because 
of the substantially stronger normal-state T dependence. 
The low-T maximum values, «;& ~ 140 and 245 W/mK 
for samples 189 and 259, respectively (and 500W/mK for 
crystal 315 by extrapolation), are the largest reported for 
anyjCuprate, exceeding the lOOW/mK of insulating Nd- 
214,0 eauLa^MW/mK observed in untwinned T c = 92K 
Y-123.yiaOt§E3 These observations indicate substan- 
tially superior and rather conventional lattice conduction 
in Y-124, confirming a similar conclusion based on pre- 
vious measurements of polycrystals.Ej 
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FIG. 1. n a (T) and K b (T) for three Y-124 crystals. Speci- 
mens 189 and 259 are lightly Al-doped. Data for Y-123 crys- 
tals (solids lines) are from Ref. |l| (T c = 92K) and Ref. || 
(T c = 60K). 



The total thermal conductivity is a sum of elec- 
tronic and lattice components, First 
let us consider the electronic contributions, n a .e and 
K&,e- The Wiedemann- Franz law, n e {T) = LqT / q-LLq = 
2.44 x 10~ 8 Wf2/K 2 ), and electrical resistivitiestJ pro- 
vide upper-bound estimates: K eja (300-RT) m 2 W/mK and 
K e . b (300K) w 6W/mK. These imply k l , q « 8W/mK and 
KL,b ~ 25 — 39W/mK, and the very large in-plane lat- 
tice anisotropy, KL,b/^L,a ^3 — 4. This is to be con- 
trasted with Y-123 where similar aiiabpes-.suggest that 
is nearly isotropic in the planes. E3c3fi3'Ej Since Y-124 
and Y-123 differ only in their CuO-chain structures, we 
must conclude that the strong lattice anisotropy in Y-124 
is directly connected with the chain-related vibrational 
spectrum. 

Anisotropy in kl can arise from either anisotropy of 
the phonon group velocities, their relaxation times, or 
both. Though optic modes can contribute to heat trans- 
port and the anisotropy, their contribution should dimin- 
ish with decreasing temperature as should any associ- 
ated anisotropyo This is contrary to our observation 
that Kb/n a in the normal-state is nearly constant or in- 
creases with decreasing T (Fig. ||) [the decrease in Kb/n a 
below T c is related to superconductivity, as discussed be- 
low] . Experimental dispersion curves arejiot available for 
Y-124, however the computed spectrumc 2 ] is quite simi- 
lar to that of Y-123, and does not imply a substantial 
in-plane anisotropy in the acoustic-mode group veloci- 
ties. These observations would suggest a large anisotropy 
in the phonon relaxation rates in Y-124, but further in- 
formation about the low-energy vibrational spectrum is 
clearly needed to address this issue. 

The rather strong influence of the light Al doping in 
the planes on k provides insight into the phonon scaU 
tering in both Y-124 and Y-123. Nuclear quadrupoleli 2 ] 
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FIG. 2. Anisotropy, K,b/^ a vs T for Y-124 crystals. Also 
shown are data for untwinned T c = 92K Y-123 from Ref. |l6| 
(solid curve) and Ref. ^ (dashed curve). 
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resonance studies of similar Y-124 crystalsreveal substan- 
tially larger 63 Cu linewidths in the Al-doped crystals, 
attributed to static disorder in the Cu(2)-apical oxygen 
bondlcngths (a typical Al-0 bondlength for octahedral 
coordination is 1.94A as compared with 2.28A for Cu- 
O). Such local structural distortions evidently represent 
a much more effective phpnon scattering mechanism than 
does the mass defecttJO introduced by the substitution. 
The local structural modifications may also explain why 
the hole concentration and T c are so dramatically sup- 
pressed - the 2% Al substitution for planar Cu reduces 
the mobile planar hole concentration by about 20%, i.e. 
a ten-fold larger reduction than would be expected if 
each Al atom filled one hole. This observation suggests 
that the suppression of mobile holes, by the combined 
effects of charge compensation, disorder-induced local- 
ization, and changes in the local electronic structure, ex- 
tends to next nearest neighbors of each Al atom in the 
planes. 

In Y-123 a similar static disorder in the apical 
bondleagths is induced by oxygen vacancies on the 
chainsc3 and there is evidence for static or dynamical 
structural distortions of the CuC>2 planes in the normal 
statc,E3 possibly related to oxygen vacancies or diffusion 
that are absent in Y-124. The implication is that these 
features give rise to substantial damping of in-plane heat- 
carrying phonons, and are the origin of the weaker In- 
dependence and in-plane anisotropy observed for k in the 
single-chain compound. 

It is of interest to compare the superconductfijig-state 
enhancement of k for Y-124 with that of Y- 123.0113 We 
define the enhancement aalJ T = — d(n s / n n )/dt\ t —>i, the 
reduced temperature derivative of the superconducting 
to normal-state thermal conductivity ratio near T c . k" is 
determined by a polynomial extrapolation of the normal- 
state data just below T c as shown in Fig. || (a). Fig- 
ure [| (b) shows the normalized data and corresponding 
slopes. 

The r values agree with those found for twinned Y- 
123 at similar p valuesJ13 indicating that for the regime 
investigated here (p < 0.11), T is largely determined by 
the change in scattering that occurs in the CuC>2 planes. 
That r a is larger for Al-free specimen 315 than for spec- 
imen 259 is consistent with the higher mobile hole con- 
centration in the planes of the former. Also of interest is 
the anisotropy: r a /Fft w 1.5 for specimen 259. Evidently 
this anisotropy is responsible for the decrease in n^/ n a at 
T <T C (Fig. |^) . A similar behavior was observed for un- 
twinned Y-123 (x ~ 0.9) [solid and dashed 
where the enhancement is also anisotropic 
a slightly smaller r a /Ffc « 1.2 — 1.5. 

Consider T in more detail. We may write 
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where i = a, b, and T e ^ and Vl% are the j-axis electronic 
and lattice slope changes, respectively. The thewnal con- 
ductivities are evaluated at T c . Theoretically^ T e is a 
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FIG. 3. (a) K a (T) and Kf,(T) showing the slope change near 
T c on an expanded scale. Specimens are represented by the 
same symbols as in Fig. 1. Solid lines are polynomial fits to 
the normal-state data, defining K n . (b) ratio of superconduct- 
ing to normal-state thermal conductivity vs reduced temper- 
ature. Solid lines are linear-least-squares fits to the data in 
the range 0.92 < t < 1, defining V = -d(n s / K n ) / dt\ t ^i. The 
data for specimen 315 (solid squares) are shifted vertically by 
0.01 for clarity. 



sum of two terms: one proportional to dA 2 /dt\ t =i and 
the other proportional to —d(T^ p /T^ p )/dt\ t= i (A and r qp 
are the real part of the superconducting gap and quasi- 
particle (qp) lifetime, respectively). In the cuprates, the 
second term predominates due to the strong suppression 
of qp scatterings at T < T c , and T e is positive. Scat- 
tering of phonons by charge carriers always produces a 
positive lattice term (Tl), but the thermal Hall conduc- 
tivity measurements of Krishana et ala indicate that the 
electronic contribution to the enhancement predominates 
in optimally-doped Y-123. Their data for twinned Y- 
123 (x ~ 0.9) yields: r a6 ~ 1.4, T e . ab ~ 12.5, and 
Ke,ab/ K ab — 9 x 10~ 2 . Thus, ^ ~ 0.3, and about 
80% of the slope change for optimally-deped Y-123 is 
electronic in origin. Recent measurementscJ suggest that 
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the slope change in underdoped Y-123 is also largely elec- 
tronic in origin, and thus it is likely that this conclusion 
also applies to Y-124. 

Given these observations, one might expect that the T 
anisotropy simply follows that of the superfluid, with the 
larger qp lifetime enhancement occurring for transport 
along the crystallographic direction characterized by the 
largest superfluid density. However, this is not the case; 
for Y-123 and Y-124 the, SJaperfluid density is 2-3 times 
larger along the chains ,Ej : EjC3 whereas F is larger along 
the planes. To get further insight into the origin of the 
r anisotropy, we examine the limiting case, T^^ = 0, for 
which Ti is determined entirely by the electronic term. 
The electronic superconducting-state slope change may 
be related to the low-frequency electrical conductivity 
via the Wiedemann-Franz law,E3 K, e i = Lia\{T , where L, 
and o\i are the Lorenz number and real part of the elec- 
trical conductivity, respectively. Substituting into Eq. |l| 
we have 
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as magnitude and sign that 
3 We first examine the see- 



Theoretically, the first term 
are highly model dependent 
ond term. For untwinned, optimally-doped Y-123, the 
microwave data of Zhang et alE2 yield, (daf a / 'g?T)|t c — 
-0.8 x 10 5 (f2 m K)" 1 and {da s lb / dT)\ Ta ~ -1.8 x lO^ft 
m K) _1 . For the normal-state data we employ da/dT = 
— (1/ p 2 )dp/dT and an average of rp{T) measurements 
on high-quality untwinned crystalsEJ: p a (T c ) ~ 75 ± 
25/^cm, dp a /dT\ Tc ~ 0.85 ± 0.25/xficm/K, p b (T c ) ~ 
38 ± 12^cm, <WdT| Tc ~ 0.35 ± 0.10/^cm/K. For the 
k data we usejli n a {T c ) = 11.5W/mK and n b {T c ) = 
16.0W/mK. With these values, the second terms in the 
expressions for T a and Ft, are, (1.2 ± 0.2lZ^j / r L(i_and 
(2.0±0.2)L b /L , respectively. ExperimentsE#jli3'PJ in- 
dicate r a = 1.5 — 1.6 and = 1.0 — 1.3. One possibility 
is that the first terms, (dLf/dT — dL™ / 'dT)\ Tc , are negli- 
gible, whence the Lorenz numbers at T c are L a / Lq ~ 1.3 
and Lb/L ~ 0.6. Though L n ^ L has been predicted 
for strongly correlated systems £3 these values for L a and 
L b are incompatible with the results of Ref. || which im- 
ply L a b ~ (0.3 — 0.5)Lo for twinned crystals. The lat- 
ter estimate is consistent -with the theoretical analysis 
of Hirschfeld and PutikkaJH which incorporates inelastic 
spin-fluctuation scattering and parameters constrained to 
fit the normal-state NMR relaxation and electrical resis- 
tivity. We thus conclude that both L a and L are less 
than Lq and that the (dLf/dT — dL™ / dT)\x c are nega- 
tive and of comparable magnitude to the second terms 
in Eq. ||. A considerable a/b anisotropy (~ 2 — 3) in 
the Li, {dLf/dT - dLf/dT)\ Te , or both is then implied. 
Calculations which incorporate inelastic as well as elastic 
scattering to yield the Li(T) near T c would be useful to 
extract further information from the T,. 
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